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Abstract: In a unified viewpoint in quantum channel estimation, we compare 
the Cramer-Rao and the mini-max approaches, which gives the Bayesian bound 
in the group covariant model. For this purpose, we introduce the local asymptotic 
mini-max bound, whose maximum is shown to be equal to the asymptotic limit 
of the mini-max bound. It is shown that the local asymptotic mini-max bound 
is strictly larger than the Cramer-Rao bound in the phase estimation case while 
the both bounds coincide when the minimum mean square error decreases with 
the order O(^). We also derive a sufficient condition for that the minimum mean 
square error decreases with the order O(^). 



1. Introduction 

In quantum information technology, it is usual to use quantum channel for send- 
ing quantum state. Since a quantum channel has noise, it is important to identify 
quantum channel. In this paper, we consider theoretical optimal performance of 
quantum channel estimation when we can apply the same unknown channel sev- 
eral times. In order to treat this problem, we employ quantum state estimation 
theory. In our setting, we can optimize our input state and our measurement [TJ 

assumed to be possible to use entanglement with reference system in the mea- 
surement process when the channel Ag with the unknown parameter 9 is applied 
n times. This setting is mathematically equivalent with the setting given in Fig. 
[5] with n = 4, which has a single input state in the large input system and a sin- 
gle measurement in the large output system. In this paper, we consider quantum 
channel estimation with the formulation given by Fig. [5] 

In the state estimation, when the number n of prepared states goes to infinity, 
the mean square error (MSE) behaves as the order O(-) as in the estimation 
of probability distribution. However, in the estimation of quantum channel, two 
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Fig. 1. Estimation scheme of quantum channel 
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Fig. 2. Simpler estimation scheme of quantum channel 



different analyses were reported concerning asymptotic behavior of MSE. As the 
first case, in the estimations of depolarizing channels and Pauli channels, the 
optimal MSE behaves as 0(^-)[IJ|3]. As the second case, in the estimation of 
unitary, the optimal MSE behaves as O ( -k ) p5U6"U7U8"U9ll 1011 1 1 j ■ In the second 
case, two different types of results were reported: One is based on the Cramer- 
Rao approach [H[S]. The other is based on the mini-max approach 6,7,8,9,10, 

The Cramer-Rao approach is based on the notion of locally unbiased estima- 
tor, and allows one to give a simple lower bound (the Cramer-Rao bound) to the 
mean square error (MSE) at a given point. The mini-max approach aims to min- 
imize the maximum of the MSE over all possible values of the parameter. The 
mini-max approach is more meaningful than the Cramer-Rao approach because 
the true value of the parameter is unknown. So, the Cramer-Rao bound is just 
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a lower bound in general, while it can be asymptotically acheived in the case of 
quantum state estimation with the n copies of the unknown state. However, the 
Cramer- Rao bound has been considered so many times in the literatures [TJ[3J|31[5] • 
The reason seems that computing the mini-max bound is a much harder problem. 
Indeed, there are only very few examples of calculations of the mini-max bound, 
and most of them are in the compact group covariant setting, in which, as was 
shown by Holevo [15] , the mini-max bound coincides with the Bayesian average 
of the MSE over the normalized invariant measure. So, many researchers [6,7,8, 
9,10,11 applied this approach to quantum channel estimation under the group 
covariance. As the most simple case for unitary estimation, phase estimation has 
been treated with the mini-max approach by several papers[9, 10, llj, and it has 
been shown that the minimum MSE that behaves as O(^). On the other hand, 
the Cramer-Rao approach suggests the noon state as the optimal input 36,37 in 
phase estimation. The later estimation scheme was experimentally demonstrated 
in the case of n = 4 [32,33 and n — 10 [34]. Also, another group [31] experi- 
mentally demonstrated an estimation protocol concerning the group covariant 
approach proposed by |35j . In phase estimation with group covariant framework, 
the asymptotic minimum MSE behaves as Of ^)[9"1[TU]. When we focus on the 
difference 9 — 9 between the true parameter 9 and the estimate 9, the limiting 
distribution concerning the random variable n(9 — 9) can be obtained through 
Fourier transform of a function with a finite domain [11) . However, the Cramer- 
Rao bound is different from the asymptotic limit of the mini-max bound. So, 
these results seem to contradict with each other. No existing research compares 
both approaches in a unified viewpoint. 

The manuscript consists of two parts. In the first part, we discuss the Cramer- 
Rao approach for channel estimation, and give a simple formula to compute a 
bound based on the right logarithmic derivative introduced by Holevo [5D]. This 
formula is based on the Choi-Jamiolkowski representation matrix [3T)ll4"0"] of the 
quantum channel, and holds under the condition that, at the given point, the 
support of the Choi-Jamiolkowski representation matrix contains the support of 
its derivative. Under this condition, we prove the additivity of the RLD Fisher 
information for quantum channels. 

The second part is about the local asymptotic mini-max bound [38 , in which 
the maximum of the MSE is taken over an interval, the width of which is sent 
to zero after computing the asymptotic limits. The obtained results are sum- 
marized as follows: (1) The local asymptotic mini-max bound is strictly larger 
than the Cramer-Rao bound in the case of phase estimation. (2) Both bounds 
coincide with each other when the asymptotic minimum MSE behaves as 0(£). 
(3) The local asymptotic mini-max bound is achievable (Proposition 3), using a 
variation of the two-step strategy [231122] . That is, there is an optimal sequence 
of estimators that achieves the local asymptotic mini-max bound at any point. 
As consequence, the maximum of the local asymptotic mini-max bound is shown 
to be equal to the asymptotic limit of the mini-max bound. 

However, the Cramer-Rao bound has a different type of achievability. That 
is, there is an optimal sequence of estimators that achieves the Cramer-Rao at a 
given point. This characteristic is a curious quantum analogue of superefficiency 
in classical statistics. In estimation of probability distribution, if we assume a 
weaker condition for our estimator, there exists an estimator that surpasses the 
Cramer-Rao bound only in measure zero points. Such an estimator is called 
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a superefficient estimator 29]. In the estimation of phase action, we point out 
that the Cramer-Rao bound can be attained only by a quantum channel version 
of a superefficient estimator that works at specific points. Since the Cramer- 
Rao approach is based on the asymptotically locally unbiased condition, we 
can conclude that the asymptotically locally unbiased condition is too weak 
for deriving the local asymptotic mini-max bound, which can be attained in 
all points. Indeed, a similar phenomenon happens in quantum state estimation 
when we use the large deviation criterion|18 . 

This paper is organized as follows. Some of obtained results are based on 
quantum state estimation with the Cramer- Rao approach. Section [2] is devoted 
to a review of the Cramer-Rao approach in quantum state estimation. In this 
section, the symmetric logarithmic derivative (SLD) Fisher information and the 
right logarithmic derivative (RLD) Fisher information are explained. However, 
the Cramer-Rao bound is obtained only by the locally unbiased condition. So, it 
is needed to discuss its relation with the estimator that works globally. In Section 
131 we treat SLD Fisher information and RLD Fisher information in the quantum 
channel estimation, and discuss the increasing order of SLD Fisher information. 
In Section 01 we give several examples where the maximum SLD Fisher infor- 
mation increases with 0(n 2 ). In Section [SJ we compare the local asymptotic 
mini-max bound and the Cramer-Rao bound. We also show the global attain- 
ability of the local asymptotic mini-max bound in the channel estimation. It is 
also shown that the Cramer-Rao bound is closely related to a quantum channel 
version of superefficiency in the phase estimation. 



2. Cramer-Rao bound in quantum state estimation 

In quantum state estimation, we estimate the true state through the quantum 
measurement under the assumption that the true state of the given quantum 
system W belongs to a certain parametric state family {pe\0 & C R d }. In the 
following, we consider the case when the number d of parameters is one. Usually, 
we assume that n quantum systems are prepared in the state pg. Hence, the total 
system is described by the tensor product space H® n , and the state of the total 
system is given by pf n . 

In this case, when we choose a suitable measurement, the MSE decreases in 
proportion to n _1 as in the estimation of probability distribution. So, we focus on 
the first order coefficient of the MSE concerning ~ . In the most general setting, 
any positive operator valued measure (POVM) M n on the total system H® n is 
allowed as an estimator when it takes values in the parameter space © C R. The 
MSE is given as 



In the quantum case, there are several quantum extensions of Fisher infor- 
mation when the state pe is diffcrentiable at 9 and (I — P) (/ — P) = 0, where 
P is the projection to the support of pg. The largest one is the right logarith- 
mic derivative (RLD) Fisher information Jg , and the smallest one is symmetric 
logarithmic derivative (SLD) Fisher information Jg . For these definitions, we 
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define the RLD Lg and the SLD L$ as the operators satisfying 
dpg _ R dpg _ 1 s s 

Then, the RLD and SLD Fisher informations are given by [T{?ll20ll2T] 

Jf := Tr P9 Lf(Lf)\ J e s := Trp e (L s 9 ) 2 . 

When the range of pg contains the range of (^-) 2 , the RLD Fisher information 
has another expression: 

When the state family {pg} is given by pg :— e z9H \u}(u\e~ l6H , the condition (JTJ 
does not hold, where H is an Hermitian matrix. In this case, the SLD Fisher 
information is calculated as follows [17]. 

4({u\H 2 \u)-{u\H\u) 2 ). (2) 

Now, we introduce the unbiased condition by 



J 8Trpf n M n (d6) = 9, M6 G 0. 



However, this condition is sometimes too restrictive in the asymptotic setting. 
So, we consider the Taylor expansion at a point 9q and focus on the first order. 
Then, we obtain the locally unbiased condition at 9q: 

J §Trpf o n M n (d6) =9 , ±j <?Trpf n M n (S)\ B=9o = 1. 

Under the locally unbiased condition at 6q, an application of Schwarz inequality 
similar to the classical case yields the quantum Cramer-Rao inequalities for both 
quantum Fisher information. 

MSE^M")^^)- 1 (3) 



MSEe (M") > —(Jg o )~ l . (4) 

Since is greater than Jg , the inequality (jU is more informative than the 
inequality ([3]). When the estimator M n is the spectral decomposition of the 

operator 9 I + + • • • + L^ (n) ), the equality in © holds, where X^ 

' »0 

is given as J®-? -1 <gi X <E) I® n ~j . Then, we obtain the following inequality 

jR > jS 

Jg S Jg . 

This inequality seems to imply that J@ is not as meaningful as Jg in the one- 
parametric case. However, as will be explained latter, provides a meaningful 
bound for MSE in the case of channel estimation. 
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In fact, in the asymptotic setting, a suitable estimator usually satisfies the 
asymptotic locally unbiased condition: 



1 



lim / 6Trpf n M n (d9) = O , lim -r n \ 9Tvpf n M n {d6) 



X) ill) 

for all points 9q. Under the above condition, using we obtain the inequality 
limsupnMSE e (M n ) > (jf)" 1 . (5) 



Further, by using the two-step method, the bound (Jg) 1 can be universally 
attained for any true parameter 9 [221125] . So, defining the Cramer- Rao bound: 

{M n } satisfies the asymptotic 
locally unbiased condition. 



C e ■= inf <^ limsupnMSE e (M") 



we obtain 



S\-l 



On the other hand, its multi-parameter case is much complicated even in 
the asymptotic setting [24, 25 ,26, 27,25]. So, this paper does not treat the multi- 
parameter case. 



3. Maximum SLD and RLD Fisher informations in quantum channel 
estimation 

In this section, we apply the Cramer-Rao approach to estimation of channel. 
In the quantum system, the channel is given by a trace preserving completely 
positive (TP-CP) map A from the set of densities on the input system H := C d 
to the set of densities on the output system JC :— C d . By using dd! linear maps Fi 

from S(H) to S(K), any TP-CP map A can be described by A(p) = Y?dx F iP F} . 
Hence, our task is to estimate the true TP-CP map under the assumption that 
the true TP-CP map belongs to a certain family of TP-CP maps {Ag}. 

In order to characterize a TP-CP map Ag, we formulate the notation con- 
cerning states on the tensor product system H®1Z, where 1Z is a system of the 
same dimensionality as ~H and is called the reference system. Using a linear map 
A from 1Z to H, we define an element \A)) of % ® 1Z as follows. 

j.k 

where {\j)ii}j=i,...,d and {\k) n}k=i,...,d are complete orthonormal systems (CONSs) 
of H and 1Z. Hence, the relation 

B®C\A)) = \BAC T )) 

holds. This notation is applied to the cases of K. (g> T-L and K.®TZ. 

Now, we focus on the matrix p[Ag] := (Ag ® id)(|/))((J|), which is called the 
Choi- Jamiolkowski representation matrix [3^130] . Then, when the input state is 
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the maximally entangled state Ittj-O) = Sj=i 77j|b') ® U)' ^ ne output state is 

2p[Ao] = (A e (gi id)(|^J))((^/|). When the matrix AA T is a density matrix 

on H, \A))((A\ is a pure state on the product system Thus, the output 

state is given as 

(A e ®id)(\A))((A\) = (I®A T )(A e ®\&){\I)){{I\){I®A) = (I ® A T )p[A e ](I <8 A). 

In the one- parameter case, we express the derivative dp ^ e ^ by £>[yle]. 

When the input state is the product state \v){v\ ® the output state is 

Ae(\v)(v\)® \u)(u\ = (I ® u ■ v T )p[Ag](I ®v ■ u t ). Since 

(?J| Tr K /o[yle]|iJ)|u)(u| = u • v t Ttk: p[A s ]v ■ w t 

= Tr K ;(/(g)M-W T )p[yl e ](/(g)lJ-W t ) = \u){u\, 

we have 

(F|(Tr K p[^])|tT) = l. 

Thus, we obtain 

Tr K p[Ag] = I. (6) 

Taking the derivative in ([S]), we obtain 

Tr K D[Ag]=0. (7) 

Now, we back to our estimation problem. In this problem, our choice is given 
by a pair of the input state p and the quantum measurement M. When we fix 
the input state, our estimation problem can be reduced to the state estimation 
with the state family {Ag(p)\9 G 0}. In the one-parameter case, we focus on the 
suprema 

J R [A e ] := sup J R [A e ,p], J s [A g ] := sup J s [A e ,p], 
p p 

where J [Ag,p] and J R [Ag,p] are the SLD and RLD Fisher informations when 
the input state is p. In particular, it is important to calculate the supremum 
J 5 [/le] which is smaller than J-^ylg]. 

When n applications of the unknown channel Ag are available, the input state 
p n and the measurement M n are given as a state on (H®lV)® n and a POVM on 
(JC®7V)^ >n . For a sequence of estimators {(p n , M™)}, we consider the asymptotic 
locally unbiased condition: 

lim f 0TvAg o ( Pn )M n {dd) = 6 , lim f 9Tv Ag(p n )M n (d§) 

n—¥oo / n— >oo qxj J 

for all points O , and denotes the MSE of (p n ,M n ) by MSEg(p n ,M n ). Assume 
that J s [Af n ] behaves as 0(n a ) when n goes to infinity. When {(p n ,M n )} sat- 
isfies the asymptotic locally unbiased condition, the inequality §5§ yields that 

limsupn" MSEg(p n , M n ) > lim sup 



J s [Af n ] 
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We define the Cramer-Rao bound: 

cm 

:= inf jumsupn a MSE 6 i(p n ,M n ) 

{p n ,M n } { n->co 



{(p n , M n )} satisfies the asymptotic 
locally unbiased condition. 



Thus, we obtain 



Since 



n 

C a [A e ] = lim sup » . (8) 



n— >oo 



J s [Af n ] 



jS [ A fn+ m ^ > jS + jS ^®m] ^ (g) 

the limit lim„_ i . 00 — n a exists. (For example, see Lemma A.l of [2I]-) Thus, 
Ci[yle] can be defined by lim^oo /S |^ n] . 

In order to treat the above values, we consider the following condition: 

(C) The range of p[Ag] contains the range of D[Ag} 2 . 

Assume that the condition (C) does not hold. When the input state is the max- 
imally entangled state l"^-0)i the RLD Fisher information diverges. So, J^Ldg] 
is infinity. 

Theorem 1 When the condition ( C) holds, 

J R [A e ] = WTncDlAeWAg^DlAo}]]. 

Proof: Assume that the input state is given by |^4))((j4| and A is an invertiblc 
matrix. Then, the range of (I <g> A T )p[Ag](I (£) A) contains the range of ((/ £§> 
A T )D[A e ]{I®A)f. 

Assume that p[Ag] is invertible. Using the formula (JXJ) , we obtain 

J R [A , \A))((A\] - Tr((7 ® A T )D[A e ]{I ® A)) 2 ((I ® A T )p[A e ](I ® A))- 1 
= Tt(I ® A T )D[A e ](I ® A)(I ® A T )D[A e ](I ® A)(I ® A)- 1 ^] -1 ^ ® ^ T ) _1 
= Tr(/ ® 3 J 4 T )£>[A 9 ]^[^ 9 ]- 1 £>[A e ] 
= TrAA T (Tr^i?[yl e ]p[yl e ]- 1 ^[yl e ]), 

where ((/ ® A T )p[yl e ](/ ® ^4)) _1 is the inverse of (I ® A T )pL4 e ](/ ® 3) on its 
range. So, the supremum of Tr AA T (Tijc D{Ag]p{Ag}~ 1 D{Ag]) with the condition 
rankX4 T = dim'H equals || Tr K D[A g ]p[A e }~ 1 D[Ag]\\. Assume that p[Ag] is non- 
invcrtible. We choose an arbitrary small real number e > 0. Similar calculations 
and the operator monotonicity of x i— > — x^ 1 yield that 

Tr((7 ® A T )D[A e ]{I ® 3)) 2 ((7 ® A T )(p[;l e ] + el)(l ® 

< TrX4 T (Tr,c #[/l e ](pL4 e ] + e/)- 1 !?^]) 

< Tr X4 T (Tr K D^]^]" 1 !)^]), (10) 
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which is a bounded value due to Condition (C). Now, we focus on two matrixes 
on the range of {I®A T )p[Ag]{I®A), {I®A T )p[Ag]{I®A) and {{I®A T )D[Ag]{I® 
A)) 2 . Since the right hand side of (fTCTI) is independent of e, the range of (7 ® 
A T )p[A g ](I ® A) contains that of ((/ ® A T )D[A S ](I ® A)) 2 . So, taking the limit 
e — > 0, we obtain 

J R [A e , \A))((A\] = Tr((7 ® A T )7}[yl e ](7 ® A)) 2 ((7 ® A T )pL4 9 ](7 ® A))" 1 
= lim Tr((7 ® A T )D[Ag](I ® ~A)) 2 ((I ® A T ){p[Ae\ + el)(l ® A)) -1 

E-5- + 

< TrA~A T (Tric D^piAg^DiAe}). 

The remaining task is to show the inequality 

J R [A e ,\A))((A\] < WTncDiAeMAor'DiAe}]] (11) 

for a non-invertible matrix A. Define (7 ® A T ) _1 and (7 ® -A) -1 as the inverses 
of 7 (8> A T and I A whose domains are the ranges of the matrixes 7 ® A T and 
I® A. Letting P and P' be the projections to the ranges of the matrixes 7 ® A T 
and 7 <g> A, we have (7 ® A T ) = P(7 ® A T ) = (7 ® A T )P' = P(7 ® A T )P', and 
(J® 3) = P'(7®2") = (7®3)P = P'(7®2")P. Then L (7® A T )p[yl e ](7(g)l) = 
P(7 ® A T )p[^ 9 ] (7 ® A)P = P(7 ® A T )PV[A e ]P'(7 ® A)P. 

In the following, we consider the case when pL^g] is invertible. The matrix 
(P' p[Ag]P')~ 1 is defined as the inverse of P'p[Ag]P' whose domain and range 
are the range of P'. Thus, 

J R [A g , \A))((A\] = TrP((7 ® A T )D[Ag](I ® 2~)) 2 P(P(7 ® A T )p[A 9 ](/ ® ^P)" 1 
= Tr(7 ® A T )P[/l e ]P'(7 ® 3)P(7 ® A)- 1 {P' p^P')- 1 (7 ® A T )- 1 P(7 ® A T )P'PL4 ] (7 ® A) 
= Tr(7 ® A T )P[/l e ]P'(PV[yl e ]P')" 1 P'P[vl e ](7 ® A) 
= Tr(P'p[^]P') _1 ^'^[^]a®3A T )P[^]P' 
= Tr(P'p[A g }P')- 1 D[A g ](I ® AA T )PL4 e ] 

<Trp[/l e ]- 1 P[^](7® AA T )P[yl e ] (12) 
= TrAA T (Tr K D^plAg^DlAg}). 

Therefore, the inequality (ITT1) holds. 

Next, we consider the case when p[A$] is non-invertible. We choose an arbi- 
trary small real number e > 0. Similar calculations and the operator monotonic- 
ity of x i-> —ir" 1 yield that 

Tr P((7 ® A T )D[A g ](I ® 2~)) 2 P(P(7 ® A T )(p[Ag] + e7)(7 ® A^P)- 1 

< TrA~A T (Trfc D[A 6 ]{p[A ] + eI)- l D[A e }) 

< TyAA t (Tt k D[A e } P [A e }- l D[A e ]), (13) 

which is a bounded value due to Condition (C). Now, we focus on two matrixes on 
the range of P, P{I®A T )p[Ag\(I®A)P and P((I®A T )D[Ag](I®A)) 2 P. Sincethe 
right hand side of (|T3"| is independent of e, the range of P(7® A T )p[Ag](I ® A)P 
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contains that of P((7 ® A T )D[A e ](I ® A)) 2 P. So, taking the limit e ->• 0, we 
obtain 

J^yle, |A»«A|] = Tr P((7 ® A T )7?[yl e ](7 ® Z)) 2 P(P(7 ® A T )p[A e ](7 ® A)P)" X 
= lim TrP((/(g)v4 T )P'[vl e ](/(g)3)) 2 P(P(/®yl T )( j o[yl e ] + e7)(7 ® ^P)" 1 

C-5- + 

< Tr AA T (Ti X DiAeMM^DiM)- 

Therefore, the inequality (fTTj) holds. ■ 

Theorem 2 When two channel families {Ag} and {Ag} satisfy the condition 
(C), then the additivity 

J R [Ag ® Ag] = J R [Ag] + J R [Ag] 

holds. 

Proof: Let K, and K, be output systems of the channels Ag and Ag. Then, the 
relation ([7]) guarantees 

Tr^ D[A e ] ® D[i fl ] = Tr/c 7JL4 e ] Tr £ £>[i fl ] = 0. (14) 
Since 7>[/hj ® i e ] = 7J[/L e ] ® p[i e ] + pL4 e ] ® 7>[i e ], Theorem Q] and dHJ) yield 

J R [Ag®Ag] 

= || Thr^^D^] (8 p[i e ] + p[A e ] ® D[yl e ])(p[yl fl ]- 1 ® p^]- 1 )^^] ® p[Aj] + pLl e ] ® £>[Aj] 

=11 ^ajeC-D^lP^l -1 - ^] ® P[M + p[A e ] ® P^lP^r 1 ^] + 27>L4 e ] ® £>[i e ])|| 
= ||(Tr JC Z^p^]- 1 ^]) ® 7 + 7 ® (Tr^ D^p^]- 1 ^])!) 
= || Tr c D^lp^]- 1 ^^]!! + || Tr £ T^p^]" 1 ^]!! (15) 
=J fl [A?] + J R [M 

where the equation (|15[) follows from the additivity property concerning matrix 
norm: 

\\x ® / + / ® y|| = |I jsc || + ||y|| 

for any two Hermitian matrixes X and y. ■ 

Corollary 1 When a channel family {Ag} satisfies the condition (C), then 

J R [Af n ] = n,J R [Ag]. 

Since nJ s [Ag] < J s [Af n ] < J R [Af n ], J s [Af n ] increases in order n under the 
assumption of Theorem[TJ i.e., J^LAf 1 ™] = 0(n). When the rank of p[Ag] is the 
maximum, i.e., ddf, this condition holds and J^/lf 1 ™] = 0(n). However, there is 
an example that does not satisfy the above condition but satisfies the condition 
(C) as follows. So, the condition (C) is weaker than the condition that p[Ag] has 
the maximum rank. 
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A channel A is called a phase damping channel when the output system K, 
equals the input system % and there exist complex numbers dk,i such that 

A {d k ,i}{p) =^2d k ,iPk,l\ k ){ l l 

k,l 

where p = J^k i Pk,l\k){l\- In this case, the state p[A] is written as the following 
form 

p[A{d k , l} ] = J2 d ^ k W( l i( l i- 

k,l 

That is, the range of p[A) is included by the space spanned by When 
a channel family {Ag} is given as a one-parameter subfamily of {A^ dk j \dk,i is strictly positive.}, 
the condition (C) holds. Therefore, there exists a channel family that satisfies 
the condition (C) but consists of non-full-rank channels. 
Further, we have the following observation. 

Corollary 2 Assume that the condition (C) holds and there exists a normal- 
ized vector u in the input system % such that W^yiq D[Ao\p[Ag\~ 1 D[Ae\\\ = 
(ulTncD^piAe^DiAe^u), [{u\p[Ag]\u),{u\D[A $ ]\u)] = 0, and [I K ®\u)(u\, p[A }} = 
0. Then, J s [Ag] = J^jle], and this bound can be attained by the input pure state 
\u)(u\ on %. That is, it can be attained without use of the reference system. 

Proof: Since [{u\p[Ag]\u) , {u\D[Ag]\u)] = 0, J s [A e , \u) (u\\ = J R [A e , \u)(u\}. So, 
it is sufficient to show that J R [A g , \u)(u\] = (u\ Tr K D[A s ]p[A e ]~ 1 D[A g ]\u). 

We consider the case when the input state is \u}(u\ €5 Remember that 

J R [A e , \u)(u\(g>\v}(v\] = J R [A g , \u)(u\}. Since [I K (g>\u)(u\, p[A g }} = 0, the equality 
in holds for A = u ■ v T . Then, 

J R [A g , \u)(u\ ® \v)(v\] = Tr (^T)( u . v Tf T[K . D [ Ae ]p[Ae]- J -D[A 9 ] 
^(ulT^DiAeUAsr'DiAsWu). 



4. Examples 

We consider the case where the condition (C) does not hold. As the simplest 
example, we consider the one-parameter unitary case, i.e., the case when A g (p) = 
e l9H pe~ t9H with an Hermitian matrix H . Using @, we obtain 

J s [A e ] = (A max (ff)-A min (ff)) 2 , 

where A max (-ff) and A m ; n (-ff) are the maximum and minimum of eigenvalues of 
H. So, we obtain 

J s [Af n ] = n 2 (A max (ff) - A min (ff)) 2 . 

/I o \ 

In particular, in the two-dimensional case, when H = I g 1 I the optimal 

input is ;^|(|0}®" + |l)® n ), which is called the noon state. This unitary esti- 
mation is called phase estimation and this estimation with the noon state was 
experimentally realized with n = 4 [321155] and n = 10[54]. 
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Next, we consider the d-dimensional system <C d spanned by {\j)}^Zo an d the 

unitary matrix X defined as X\j) := \ j + 1) mod d. Using a distribution {pj}fco 

and a real diagonal matrix H with diagonal elements {h^Z^, we define the 
TP-CP map A 8 by 

MP) :=Y.V^^ m pe- im X-K 

3=0 

This TP-CP map can be regarded as the stochasitc application of the unitary 
X 3 after the application of the unitary e l8H . 

Let h a and hi, be the maximum and the minimum eigenvalues of H. Using two- 
dimensional reference system spanned by \0) R and \1}r, we choose the following 
input state: 

|* n ) :=-L(|a)®»|0) a +|&)®»|l>*). 

In this case, as the first step, we apply the following measurement {Mj}: 

Mj :=\j + (a,...,a))(j + (a,...,a)\®\0) R R (0\ 

+ \j + (b,...,b))(j + (b,...,b)\®\l) R R (l\. 

When the outcome of this measurement is j , the resulting state is the pure state 

^(e inha \j + (a,.. .,a))\Q) R + e inh »\j + (b,.. .,b))\l) R ). (16) 

The SLD Fisher information of the above family is n 2 (h a — h ) 2 . Therefore, 
since the maximum SLD Fisher information behaves as 0(n 2 ) at most, J s [Af n ] 
behaves as 0(n 2 ). 



5. Local asymptotic mini-max bound 

In this section, we consider the relation between the discussion in the previous 
section and estimating protocols in a different viewpoint. Consider the phase 
estimation with inputing the noon state ^(|0}® n + 11)®"). Then, the output 

state is -^(e™ 6 / 2 ^)®" + e _<Tie / 2 |l)® n ). In this case, we cannot distinguish the 

parameters and 6+ For example, when we apply measurement {-^(|0)®"± 

\l)® n )}, the probability with the outcome -L(|0)® n + ll)®") equals cos 2 n9/2, as 

is shown in Fig[3J and the Fisher information equals n 2 . Even if the parameter 
9 is assumed to belong to (0, Tv/n], we cannot distinguish the two parameters 
9 = 7r/3n and 9 = 2ir/3n with so high probability because we have only two 
outcomes. In order to distinguish two parameters 9 — ir/3n and 9 = 27r/3n in 
this measurement, we need to repeat this measurement with several times, e.g., 
k. Since the number of application of the unknown unitary is N := kn, the error 
behaves as r^, which is different from O(j^). So, we cannot conclude that the 
above method attains the order O(-jb) concerning MSE. Therefore, we need to 
discuss what a bound can be attained globally, more carefully. 
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Fig. 3. Phase estimation with noon state n = 20 



For this purpose, we focus on an e-neighborhood Ug, c of and define the local 
asymptotic mini-max risk [58]: 

C ol [A 0o ,{(p n ,M n )}] := limlimsup sup n a MSE# (p n , M n ) 

and the local asymptotic mini-max bound: 

C a [A 8o ] := inf C a [Ag ,{{p n , M")}j. 

{(p„,M»)} 

Concerning the local asymptotic mini-max bound, we have the following two 
propositions. 

Proposition 1 When C a [Ag n ] is continuous and the convergence is com- 
pactly uniform, 

C a [Ag ] > C a [Ag ]. (17) 

Proof: For any 5 > 0, we choose an integer TV and an e-neighborhood Ug 0j€ 
satisfying 

C a [Ag ] +S> n a MSV e (p n ,M n ), Vn > N, V0 e U 8o , e . (18) 
We introduce two quantities 



r, n (6) := / 0TrAg(p n )M n (d6) 

v n (6) := [(§-j] n (e)) 2 TrAe(p n )M n (d§). 
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Then, we obtain 

MSE e (p„,M") = v n {6) + ( Vn (9) - 9) 2 . 

Let Jg^ n be the Fisher information of the distribution family {Tr Ag(p n )M n (d6)\6 € 
<9}. This quantity is smaller than J 5 'L4f m ]. Deforming the classical Cramer- Rao 
inequality, we obtain 

v n {6) > { -^f—L. 
Je,n 

Thus, we obtain 

MSE e (p„, M n ) > V f ' + (Vn(9) - 9) 2 . (19) 

J0, n 

As is shown later, for any 6 > 0, there exists a sufficiently large integer N 
satisfying the following. For any n > N, there exists 9 n € Ug . e such that 

dT]M >l-S. (20) 
d9 

Using (fTUj) . we obtain 

MSE e „ {p n , M n ) > 1 f ' . (21) 
Je n ,n 

Take the limit n — > oo. Then, the continuity of C a L4g ], the compact uniformity 
of the convergence J8]), (fl8)l . and (|2ip imply that 

C a [A eo } + 5> (1 - 8fC a [Ae }. 

Taking the limit 6 — > 0, we obtain (fTT)) . 

Finally, we show the existence of the integer N satisfying the above condition 
given by ((20)) by using reduction to absurdity. Assume that for any S > 0, there 



exists a subsequence such that dVn ^g 6 ^ < 1 — 5 for any 9 6 e - Thus, 

[ Vnk (9 + e/2) -{9 + e/2)] - (0 - e/2) - (0 - e/2)] 
=?7 „ fe (9 + e/2) - r?„ fc (9 - e/2) - e < -eS. 

Then, 

max{|r?„ fc (0 + e/2) -(9 + e/2)|, 1^(0 - e/2) -(9- e/2)|} > y . 
That is, 

max{(r K (0 + e/2) -(9 + e/2)) 2 , ( Vrik (9 - e/2) - (9 e/2)) 2 } > ff. 
Using ((T9J), we obtain 

2 ?2 

Since MSEg +e / 2 (Pn, M n ) behaves as O(^), we obtain contradiction. 
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Proposition 2 Assume that E :— sup^g, \6\ < oo. When the order parameter 
a equals 1 and C a [Ag ] is continuous, 

C 1 [A eo ]=C l [A 0o ] = lim 



J S \A\ 



Proof: As is mentioned in Section[3J © guarantees the convergence of lim n _ s . 00 , \, . 

It is enough to show the inequality CiL4e ] > Ci[-^e ]- 

For an arbitrary real number 8 > and an arbitrary integer to, let {p m , M m } 
be a locally unbiased estimator at do such that 

' +S>MSE e<) (p m ,M m ). (22) 



JS[A% 

We define another coordinate n(0) by 

v{6) := / §Tr A e (p n )M n (d§), 



and denote the MSE concerning the parameter n of an estimator (p n ,M n ) by 
M$E rt rg\(p n ,M n ). For any 5 > 0, we choose an integer to, a sufficiently small 
e'-neighborhood Ug 0:<L such that 



Cx[A 6o ] ^ 1 



- 6-6' - K ' 

for V<9,6>' € C/ e0i£ . 

Next, let (p' nm , M nml ) be the estimator given as the average value of n times 
applications of the estimator (p m ,M m ) concerning the original parameter 8. 
Then, we can choose a sufficiently large number n satisfying the following: When 
the true parameter is 9 , the estimate of (p' nm , M nm ') belongs to Ug . t with the 
probability higher than 1 — 8. The second inequality in guarantees that 

-MSVg(p m ,M m ) + 8E 

n 

= MSE 9 (p' nm , M nml ) +8E> -^-^ MSE, (e) (^ m , M nm '). (25) 



Thus, ([22 ]) . (|23 |) . and (|25 ]1 imply that 

(1 + sfi^M + 2 A +5E )> MSV v{eo) (p> nm ,M™'). 
nm n 

Since MSE v( g a) (p' nm ,M nm ') is continuous and (I + 8f{^^ + S(E + 2)) > 

(1 + 8) 2 ( Cl ^ ^ + 77 + 8E), we can choose a sufficiently small number < e' < e 
such that 

(1 + 8) 2 {Q^A + 6 (E + 2)) > MSE v(e) (p' nm , M nm ') (26) 
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for 9 £ U 6o , e >. 

Let (pi nm ,M lnm ) be the estimator given as the average value of I times 
applications of the estimator (p' nm , M nm ') concerning 77. Since the estimator 
(p' nm , M nm ') is unbiased concerning the parameter 77, 

MSE m (p? nm> M<™") - MSE " w(P ;""' M " m,) . (27) 

We choose a sufficiently large number I satisfying the following: When the true 
parameter is 9 € Ug 0<e i / 2 , the estimate 77 of (p[' nm , M lnm ) belongs to U$ 0l£ ' with 
the probability 1—pi, where the probability pi exponentially goes to as I — > 00. 
The first inequality in ((24)) guarantees that 

MSE„ (e) (d' nm , M lnm ") + E Pl > MSMPi nm ,M lnm "), V0 G U 0( ^ /2 . 



{1-5) 



(28) 



Thus, the relations (|26 ]> . f|27 |t . and l|28 j) imply 

^^(^^+5(E + 2)) + E Pl >MSEe(pl m ,M lnm ''), WeC/ v/2 . 



(29) 

Taking the limit I — > 00, we obtain 

i i \ + S x ) .l (C 1 [A eo }+nm6(E + 2)) > lim Inm MSE^'™, M lnm "), V0 € U 0o . e , /2 . 

Finally, using the above n and to, we define a sequence of estimators {(p' fc ",Af fc )} 
by the following way. For given fc, we choose maximum / such that Inm < k. 
Then, the estimator {p'^M k "') defined as {p'( nm , M lnm "). In this definition, we 
only use Inm applications, and the remaining k — lnm applications are discarded. 
So, we obtain 

S! _ A?2 [^0] + nmS[E + 2)) > lim sup k MSE 9 (p' fc ", M k "'), V6 € U eo>€ , /2 , 
which implies that 



(i-sy 



(C 1 [A eo ]+nmS(E + 2)) > Ci[A 9o , {(/>£', )}£° =1 ] > Ci[A 9o ]. 



Since (5 > is arbitrary, CiLigJ > CiL4g ]. ■ 
Now, remember that the Cramer- Rao bound can be attained by using the two- 
step method in the case of state estimation. By using the two-step method[22, 
123] . the local asymptotic mini- max bound C Q [/l0 o ] can be attained at all points 
9 as follows. 
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Proposition 3 Assume that E :— sup eg0 \8\ < oo and Ci[Ag ] is continuous. 
For any 5 > 0, there exists a sequence of estimators {(p n , M n )} such that 

C a [A e ,{(p n ,M n )}} < C a [A e ] + 5 (30) 

for all points 9. Further, when the parameter space is compact, 

limn" min maxMSEfl(p„, M n ) = maxC„ \A e }. (31) 
n-Kx> (p„,m«) see see 

Proof: We use a two-step method slightly different from [25] . Before apply- 
ing the unknown channel Ag, for any real number 5 > 0, we choose an ej- 
ncighborhood C/g i;Ci satisfying the following three conditions: (1) UiUg i , tl = 0. 
(2) For any there exists a sequence of estimators {(/9 n (0.j), M n (0j))} such 
th&tC a [Ag t }+5/2 > lim^olimsup^^ n a sup e&Ue _^ MSE e (p„(^),M™(^)). (3) 
sup, :(V , . CUA?,] < C Q L4 e ] +<J/2. 

We divide n applications of the unknown channel Ag to two groups: The 
first group consists of y/n applications and the second group consists of n — y/n 
applications. In the first step, we apply a POVM M to the first group. This 
POVM M is a POVM on the single system H and is required to satisfy that 
Jg is non-degenerate at all points 9. Based on y/n obtained data, we estimate 
which e^-neighborhood Ug ir€i contains the true parameter, and obtain the first 
step estimate 9~ v The error probability Pg tTl of this step goes to exponentially, 
i.e., Pg >n behaves as e~ Cv/ ™, where c depends on 9. 

In the second step, we apply the estimator (p n _^(9i),M n ~^{9i)) to the 
second group, and obtain our final estimate from the outcome of the estimator 
(p n -y/n(Oi),M n -^ K (0i)). We express this estimator by (p n ,M n ). Its MSE is 
evaluated as 

MSEg( Pn , M n ) < EP e . n + (1 - P e>n ) sup MSEflOw^fli), M n -^ I (9 l )). 

v.eeu a . H 

(32) 

Since n a EPg, n goes to 0, we obtain 

limsupn Q MSE e (p„,M n ) < sup C a [A 8i ] + 5/2 < C a [Ag] + 5. 

Thus, we obtain (f3T)l) . Further, the relation (|3"2"j) yields that 

supMSE e (p„,M") 
eee 

<supEPg. n + sup(l~Pg. n ) sup MSEg(p n _^(9 t ),M n -^(9 l )). 

The compactness of guarantees that sup ee6 > n a EPg tJl —> 0. Thus, 

limsupn" sup MSE e (p„, M n ) < sup sup C a [A Bi ] + 5/2 < supC a [A e ] + 5. 

n-too gee 8e&i-9eu SitH eee 



Since the part > of (|3"Tj) is trivial, we obtain (pri) . 
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Proposition [3J holds even when we replace the MSE by a general error func- 
tion R(9,9) for one-parametric family satisfying the following conditions: (1) 
the relation R(9, 6) = (6 — 9) 2 holds with a local coordinate when 9 is close 
to 9. (2) the maximum of R(9, 9) exists. Therefore, we can apply Proposition 
to the following case: Assume that the one-parameter channel family {Ag} 
has a compact group covariant structure, that is, its parameter space is given 
as an interval [a, b) and there is a unitary representation Ug of M such that 
Ue> Ag(p)Ug, = Ag+0'(p). The error is given by mmk e z{9 + k(b — a) — 9) 2 instead 
of the square error (9 — 9) 2 . In this case, due to the group covariance, C Q [yle] 
does not depend on the true parameter 9. Application of Proposition [3] implies 
that the global min-max error behaves C a [Ag]^. 

In the phase estimation case, the unknown parameter 9 belongs to [0,27r), 
and the minimum of the worst value of the average error maxg MSEg(/9„, M n ) 
behaves as ^ [ MTUIfTT] . That is, the leading decreasing order is 0(1/ n 2 ) and 
the leading decreasing coefficient is ir 2 when we apply the optimal estimator. 
Proposition [3] implies that maxg C2[Ag) — ir 2 . Since C2[/le] does not depend on 
9 due to the homogenous structure, we can conclude that C^vle] = tt 2 - So, the 
equation J s [Af n ] = n 2 implies the equation C^Lde] = 1. Hence, the Cramer- 
Rao bound C^^lg] cannot be attained globally in this model. However, it can be 
attained in a specific point in the following sense. 

Proposition 4 Assume that E := sup egg) \9\ < oo and C\[Ag Q ] is continu- 
ous. For any S > and any 9q G O, there exists a sequence of estimators 
{(Pn,e , Mg )} satisfying the asymptotically locally unbiased condition and the 
relations: 

lim sup n a MSE e [p nfio , M e ™ ) < C a [Ag] + 8, ^ 9 Q 
lim sup n a MSE 9o (p n .g , M£ ) < C a [Ag ] + S. 

n— >oo 

In estimation of probability distribution, there exists a superefficient estima- 
tor that has smaller error at a discrete set than the Cramer- Rao bound [25]. Since 
such a superefficient estimator cannot be useful, many statisticians think that 
it is better to impose a condition for our estimators for removing superefficient 
estimators. In this classical case, if we assume the asymptotic locally unbiased 
condition, we have no superefficient estimator. Proposition 2] means that even if 
the asymptotic locally unbiased condition is assumed, there exists an estimator 
that behaves in the similar way to a superefficient estimator in the case of unitary 
estimation. So, we call such an estimator a q-channel-superefficient estimator. 
That is, a sequence of estimators {(p n , M n )} is called q- channel- superefficient 
at 9 with the order ^ when limsup,^^ n a MSEg(p n , M n ) < C a [Ag]. Hence, in 
order to remove the q-channel-superefficiency problem, it is better to adopt the 
bound C a [Ag] as the criterion instead of C a [Ag] . 

Proof: We choose ^-neighborhoods Ug uCi in the same way, and define the neigh- 
borhood Ug a i . We apply the same first step as PropositionUJto neighborhoods 

{Ug iiCi }i U {Ug a> i } ; and obtain the first step estimate 6*-. When the first step 

estimate 9^ is not 9q, we apply the same method as Proposition [3] in the second 
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step. When the hrst step estimate #j is 9q, we apply the asymptotically locally 
unbiased estimator whose MSE behaves as (C^-Ag] + 5)/n a asymptotically. ■ 
Further, since there exists an asymptotically locally unbiased estimator that 
surpasses the bound C a [Ag] , the asymptotically locally unbiased condition is too 
weak for deriving the local asymptotic mini-max bound, which is more mean- 
ingful. In order to avoid this problem, it is sufficient to impose the following 
condition: 

(CU) The limit lim„_ i . 00 n a MSE$(p n , M n ) exists for all 8 and this convergence 
is compactly uniform concerning 9. 

Under the condition (CU), Iim„_>oo n a MSEg(p„, M n ) is continuous concerning 
6, and 

limn" sup MSE e (p n ,M n )= sup lim n a MSE e (p n , M n ). 
n->oo eef / fl0ie eeu eo , e n ^°° 

Thus, 

lim n a MSEe ( Pn ,M n ) = lim sup lim n a MSE e (p n ,M n ) 

n— >-oo Q^JJq n— >-oo 

= lim limn" sup MSEg(p n ,M n ) = C a [Ag ,{(p n ,M n )}\. 
Therefore, we obtain the following corollary. 

Corollary 3 When a sequence of estimators {(p n , M n )} satisfies the condition 
(CU), 

lim n a MSE e (p n ,M n ) > C a [Ag}. 

n— > oo 

Therefore, the condition (CU) is better in estimation of quantum channel than 
the asymptotically locally unbiased condition. 

Finally, we consider the relation with the adaptive method proposed by 
Nagaoka|16 . In this method, we apply our POVM to each single system H, 
and we decide the fc-th POVM based on the knowledge of previous k — 1 out- 
comes. In this case, Fujiwara j3U] analyzed the asymptotic behavior of the MSE 
of this estimator. Now, we consider the case of nm applications of the unknown 
channel Ag. In this case, we divide nm applications into n groups consisting of m 
applications. When we apply the adaptive method mentioned in Fuiiwara|30j to 
these groups, the MSE of this estimator behaves as — s 1 ^ , which is close to 

C °i^ . So, when a > 1, this method cannot realize the optimal order 0{ * \ a ). 



6. Discussion 

We have compared the Cramer- Rao bound C^Jyle] and the local asymptotic 
mini-max bound C a [ylg] in quantum channel estimation, which contains quan- 
tum state estimation. When the model has group covariant structure, the local 
asymptotic mini-max bound C a [Ag] coincides with the limit of the global mini- 
max bound. We have also shown that both bounds C^yle] and Ccfylg] coincide 
in quantum channel estimation when the maximum of SLD Fisher information 
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J s [Af n ] behaves as 0(n). The case of state estimation can be regarded as a 
special case of this case. That is, the conventional state estimation has no dif- 
ference between both bounds. However, we have shown that the Cramer- Rao 
bound CcJ-Afl] is different from the local asymptotic mini-max bound C^Lle] in 
the phase estimation. So, we can conclude that the local asymptotic mini-max 
bound C Q L4g] is more meaningful and does not necessarily coincide with the 
Cramer- Rao bound C a [Ag] . 

In order to clarify the asymptotic leading order of J s [vlf 1 ™] , we have derived 
the condition (C) as a sufficient condition for J s [yl|" 1 ] = 0(n). That is, the 
condition || Tr^ DlAg^Ag]' 1 D[A S ] \\ = oo is a necessary condition for square 
speedup. This condition has been derived from the following two facts. One is 
the supremum of the RLD Fisher information satisfies the additive property. 
The other is the RLD Fisher information is an upper bound of the SLD Fisher 
information. This, the supremum of the RLD Fisher information is the upper 
bound of the regularized supremum of the SLD Fisher information, which equals 
the inverse of the Cramer-Rao bound. However, it is an open problem to clar- 
ify whether this upper bound can be attained by the regularized SLD Fisher 
information. 

Further, Fujiwara and Imai [2] and Matsumoto [14] also obtained another 
sufficient condition. Since the relation with their conditions is not clear, its clar- 
ification is an open problem. Our condition (C) trivially contains the case when 
the state p[Ag] is a full rank state on the tensor product system while it is not so 
easy to derive the above full rank condition from Fujiwara and Imai's condition. 
Further, we have also obtained another example for J^^g 1 ™] = 0(n 2 ) under the 
condition (C). This example is a larger class than the unitary model. So, we can 
expect that J s [Af n ] behaves as 0(n 2 ) if the condition (C) does not hold. This 
is a challenging open problem. 
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Lemmas needed for Theorem [l] 

Lemma 1 Any strictly positive definite matrix A and any projection P satisfy 
the inequality 



A- 1 > (PAP)- 1 , 



(33) 



where (PAP) 



-l 



is the inverse matrix with the domain P. 
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Proof: Let R be the operator norm of the matrix PA(I-P). Lemma[3]guarantees 
the inequality 



i? 2 

PA(J - P) + (I - P)AP < eP + — (I-P) (34) 

e 



for any e > 0. Thus, 

A < PAP + (I- P)A(I -P) + eP+ — (I - P) 

e 

=P(A + eI)P + (I - P)(A + - P) 

Since the function x i— > — x _1 is operator monotone, 

A' 1 > (P(A + e^P)- 1 + {{I -P)(A+ ^-/)(J - P))- 1 > (P(A + e7)P)" 1 . 

Taking the limit e — > 0, we obtain ■ 

Lemma 2 Let A be a positive semi-definite matrix and P be a projection. Then, 
the inequality 

R 2 

PA(I - P) + (I - P)AP < eP + — (I - P) (35) 

e 

holds for any e > 0, where R is the operator norm of the matrix PA{I — P). 
Proof: Choose an arbitrary normalized vector u. Let t be ||Pu|| 2 . Then, 
(u\PA(I - P) + (I - P)AP\u) < 2ViVl~tR 
<te+(l-t)^ = (u\eP+^{I-P)\u), 
which implies (|35p . ■ 
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